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Summary 

A versatile large-signal, two-dimensional computer program 
is used by NASA to model coupled-cavity traveling-wave tubes 
(TWT’s). In this model, the electron beam is divided into a 
series of disks, each of which is further divided into axially 
symmetric rings which can expand and contract. The trajec- 
tories of the electron rings and the radiofrequency (RF) fields 
are determined from the calculated axial and radial space- 
charge, RF, and magnetic forces as the rings pass through a 
sequence of cavities. By varying electrical and geometric 
properties of individual cavities, the model is capable of 
simulating severs, velocity tapers, and voltage jumps. The 
calculated electron ring trajectories can be used in designing 
magnetic focusing and multidepressed collectors. The details 
of using the program are presented, and results are compared 
with experimental data. 

Introduction 

A one-dimensional computer model for a coupled-cavity 
traveling-wave tube was developed at NASA in the 1970’s 
(refs. 1 and 2) and was used to design the traveling-wave tube 
for the NASA Communications Technology Satellite (CTS). 
In 19S4, O’Malley (ref. 3) revised the model by simulating 
the electron beam with axially symmetric rings of electric 
charge which could expand or contract in the presence of a 
magnetic field. In the present study, the model has been further 
revised to allow (1) a tube with up to 200 cavities to be 
simulated, (2) magnetic pole pieces to be of different radii than 
the tunnel tips, and (3) magnetic periods to be independent 
of cavity periods in periodic permanent magnet (PPM) 
focusing designs. O’Malley referred to the model as three- 
dimensional since the electron rings travel in the axial 
direction, expand and contract in the radial direction, and rotate 
in the azimuthal direction. However, since the rings are axially 
symmetric and there is no azimuthal dependence, it is more 
accurate to refer to the model as two-dimensional. 

In this model, the electron beam is divided into a series of 
disks, each of which is further divided into axially symmetric 
rings. The trajectories of the electron rings and the RF fields 
are determined from the calculated axial and radial space- 
charge, RF, and magnetic forces as the rings pass through a 
sequence of cavities. Each cavity has individually entered 
geometric and electrical parameters. This enables the model 


to simulate severs, voltage jumps, and velocity taper designs. 

With the original one-dimensional model, Connolly and 
O’Malley (ref. 1) showed that the calculated small-signal gain 
as a function of frequency for the CTS TWT compared well 
with the experimental data. Even better agreement is expected 
with the two-dimensional model because the radial space- 
charge, RF, and magnetic forces can be simulated. Further- 
more, modeling both the axial and the radial component of 
the electron beam trajectory is essential in designing magnetic 
focusing or a multidepressed collector at the tube output. 

In this report, a description of the physical formulation of 
the model is reviewed, and the mechanics of utilizing the 
computer program are detailed and updated. The accuracy of 
the two-dimensional model is verified by comparing its results 
with those of the one-dimensional model and with the 
experimental data for the CTS coupled-cavity traveling- wave 
tube. 

Symbols 

a tunnel radius 

a inner radius of magnet stack 

B magnetic field 

B 0 axial magnetic field at pole piece radius and midpoint 
of gap 

b o initial beam radius 

c velocity of light 

E electric field 

/ frequency 

g half-length of K th magnetic gap 

7 0 beam current 

j ring index (j = 1, 2, 3 refers to inner, middle, and 

outer ring) 

L cavity length 

i K half-length of cavity gap 
m ring mass 

m electric field shape parameter 

N c number of cavities in TWT 

N d number of disks in a beam wavelength 
N z number of axial divisions in cavity 


i 



p 

magnetic period 

Pf,K 

average power flow in forward direction 

Pin 

input power 

q 

ring electric charge 

R 

number of rings per disk 

r 

centroid radius 

r 

ring radial velocity 

r R, 0 

initial radius of outermost ring 

Hn J 

inner radius of 7 th ring 

*out J 

outer radius of / h ring 

ti(Zn,k) 

time of arrival of I th ring at axial position z n j 

u 

relativistic factor 

u 

beam velocity 

V b ,K 

backward propagating voltage in £ th cavity 

Vk 

amplitude of voltage across gap in K 111 cavity 

V f.K 

forward propagating voltage in K cavity 

y< o 

beam voltage 

A V biK 

backward induced voltage in K th cavity 

A V fX 

forward induced voltage in cavity 

V P 

wave phase velocity 

%k 

total interaction impedance for cavity 

z 

axial distance 

z 

ring axial velocity 

Az* 

length of axial division in AT 111 cavity, L K /N Z 

a +j&\ 

propagation constant 

\ 

beam wavelength 


magnetic field shape parameter 

<P 

ring azimuthal velocity 

^(r 9 z) 

axial magnetic flux through ring 


cathode flux at outermost ring 

CO 

angular frequency 

Subscripts 

b 

backward 

f 

forward 

in 

inner 

K,m 

cavity numbers 

out 

outer 

RF 

radiofrequency 

r 

radial 


SC space charge 

VJ voltage jump 

z axial 

<p azimuthal 

0 initial 


Description of Model 

The theory and mathematical derivation of the governing 
equations for the two-dimensional model is documented in 
detail by O’Malley (ref. 3). In this section, the procedure for 
the solution of these equations is reviewed and summarized. 

Figure 1 illustrates the model. The electron beam is divided 
into a series of disks, each of which is divided into a maximum 
of four axially symmetric rings. The beam trajectory is 
described by following the rings contained in a single beam 
wavelength, which is equal to u 0 /f, where u 0 is the initial 
beam velocity and /is the frequency of the RF signal. The 
axial thickness of all rings remains equal and constant at 
\/N d , where \ € is the beam wavelength and N d is the number 
of disks in a beam wavelength. The rings expand and contract 
radially and may penetrate each other in both the radial and 
axial directions. When the rings enter the tube they have the 
same charge and the same cross-sectional area. The axial and 
radial thickness of a ring is taken into account only in the 
modeling of beam interception and in the calculation of space- 
charge forces. For all other forces, it is assumed that the charge 
and mass of a ring are concentrated at the ring’s centroid 
radius. The initial centroid radii are given by 



where R is the number of rings per disk and b 0 is the initial 
beam radius. The inner and outer radii of the rings are given by 
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Figure 1.— Model of coupled-cavity traveling-wave tube. 
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Thus 


(3) 


ti(Zn,K) — 4 ^ ^ — ^l,x) 


* 


n=l,...,N z (4) 


Note that the inside “ring” is actually a disk; that is, its inner 
radius is zero. 

The tube body is modeled as a conducting tunnel of radius 
a divided axially into a series of discrete cavities, where the 
length of the K * cavity is L K . In the center of each cavity is 
a gap of length 2^. Impressed across the K th gap is a 
complex voltage V K ^\ 

The procedure for calculating the electron beam trajectories 
and cavity voltages in the model is summarized in the following 
steps: 

(1) The input data for the beam and cavities are read into 
the model. The details of entering the input parameters are 
presented in the next section. 

(2) Space-charge tables are calculated. The space-charge 
force on a ring of centroid radius r x and electric charge q x 
due to a ring of centroid radius r 2 and electric charge q 2 
located an axial distance z away is given by equations (39) 
to (41) of reference 3. Using these equations for calculating 
space-charge forces on each integration step would be very 
time consuming. Thus space-charge tables are calculated for 
discrete values of r u r 2 , z, and j ( j is the ring index; for 
example, if there are three rings per disk, y =1,2, and 3 refers 
to an inner, middle, and outer ring, respectively). During the 
computer simulation, the space-charge force on one ring due 
to another is calculated by using linear interpolation with the 
appropriate table. 

(3) For the cavity, tables of coefficients are calculated 
for the RF electric field, the voltage-jump electric field, and 
the beam- wave interaction field. These coefficients are 
evaluated for an array of r and z values and stored in tables 
when cavity K is entered. The tables are interpolated in a later 
step to calculate the forces. 

In reference 3, the RF electric field coefficients are given 
by equations (59), (62), and (63); the voltage-jump coefficients 
by equations (74), (76), and (77); and the beam- wave 
interaction coefficients by equations (59) and (121). 

(4) A first approximation for the induced voltage in cavity 
K due to the electron beam is calculated. The total induced 
voltage in the A^ 1 cavity consists of a forward induced voltage 
component A V f K , which propagates in the forward direction, 
and a backward induced voltage component AV b Kt which 
propagates in the backward direction. A first approximation 
to AV f)K is obtained by assuming that the rings have constant 
velocities in the cavity. 


where (z„,at) is the time of arrival of the / th ring at the axial 
position z n> K and dt^zx^/dz = Vz(z\ tK ) is the reciprocal of 
the velocity of the 1 th ring at the axial position Z\ tK . The axial 
position is defined by 


Z>n,K — 


K- 1 

D L m + 

m = 0 



i^=l,...,^ 


(5) 


where N z is the number of axial divisions in the cavity, N c 
is the number of cavities in the TWT, L m is the length of the 
m th cavity, and A zk is the length of an axial division in the 
K th cavity, where Az* = L K /N Z . With the approximation of 
equation (4) the forward induced voltage AV f K is calculated 
from equations (126) and (127) of reference 3. 

For a forward propagating wave in the passband of a uniform 
structure in the absence of the electron beam, Floquet’s 
Theorem gives 


V f ,K = e- (a+ ^ ) % K . i ( 6 ) 


where V f K is the forward propagating voltage in the 
cavity, a + j$\ is the complex propagation constant, and L 
is the cavity length. When the voltage induced by the beam 
is added and the differences in impedances between cavities 
are taken, the forward propagating voltage is 


V f,K = 



Vf K _xe~^ aL ^f K ~ 


+ A V ftK (7) 


where the factor ^JZ K /Z K _ X ensures that power flow is 
conserved. In the first cavity K = 1, 


Vf. 1 = + A^/,1 


( 8 ) 
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where Z x is the total interaction impedance and P m is the input 
power. 

(5) In this step, the ring trajectories are calculated from the 
equations of motion. 

The axial and radial equations of motion are respectively 



where q is the ring electric charge; m is the ring mass; 


0 = 



z 2 + r 2 + r 2 <p 2 


(ID 


where c is the velocity of light; 


where c is the velocity of light, \q/m\ is the charge to mass 
ratio, and V 0 is the beam voltage. The initial azimuthal 
velocity of each ring is 



2tt4o 


(15) 


where \p c is the cathode flux at the outermost ring and r R 0 is 
the initial radius of the outermost ring. 

The cavity is divided into N z equal parts of length A z K . By 
using finite differences with z as the independent variable, 
equations (9) and (10) are integrated to the end of the cavity. 
At the end of each integration step, the solution is checked 
for beam interception. If the outer radius of a ring exceeds 
the tunnel radius a , then the outer radius is redefined to be 
a , and the portion of the ring located beyond r — a is assumed 
to be intercepted. The centroid radius for the remaining partial 
ring is redefined by 



(16) 



\p(r,z) - ^(/q,Zq) 
2tt r 2 


( 12 ) 


where \ l/(r,z) is the axial magnetic flux through the ring; B z , 
B n and B ^ are axial, radial, and azimuthal components of the 
magnetic field, respectively; £rf , z and E RF r are axial and 
radial RF electric field; z and E w J>r are axial and radial 
voltage-jump electric field; E SC z and E SCr are axial and 
radial space-charge electric field; and U z and U r are axial and 
radial relativistic factors. 

The initial axial velocity for each ring in cavity K = 1 is 



where b 0 is the initial beam radius, and the initial beam 
velocity is 


(6) After the beam ring trajectories have been calculated 
for a cavity, a better approximation to the induced voltage is 
obtained from equations (126) and (127) of reference 3. If the 
new value of Vf K differs from the old value by more than a 
predetermined amount (expressed by TOLDV in input), step 
5 is repeated until the value converges. Then the backward 
induced voltage is calculated from equations (124) and (128) 
of reference 3. The AV b K 's are stored for use in the next pass 
through the tube. 

(7) Steps 3 through 6 are repeated for each cavity in 
succession to the end of the tube. If cavity K + 1 has the same 
input parameters as cavity K , the tables in step 3 do not need 
to be recalculated. 

(8) After AV b K for the last cavity K = K 2 is calculated, the 
backward voltages are determined for all the cavities from the 
stored AV b ^’s. For the last cavity 


V b,Ki = AK, 


b,K 2 


(17) 



(14) 


Then proceeding backwards to the first cavity K — 1 
V b .K - A V bil 


^K+ 1 


K = K 2 - 1, K 2 — 2,...,1 (18) 
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(9) After all the V b K J s have been calculated, a second pass 
is made through the tube. The second pass yields a new set 
of AV b >K that can be used for calculating a set of V b K for a 
third pass. As many passes as necessary to obtain convergence 
are performed. 

Model Input 

Three types of magnetic focusing fields can be modeled: 
(1) uniform solenoid focusing, (2) single-period periodic 
permanent magnetic (PPM) focusing, and (3) double-period 
PPM focusing. 

For uniform solenoid focusing, the magnetic field is given by 


B z (r,z) = B 0 

(19) 

B r (r,z) = 0 

(20) 


magnetic shaping factor is Ji (MU). The axial magnetic field 
at the pole piece radius and midpoint of the gap is B 0 (BO). 
The magnetic field components B z (r,z) and B r (r,z) are then 
given by 

oo 

B z (r,z) = £ a n I 0 (k n r)cos(k„z) (21) 

n= 1 


B r (r,z ) = £ a„/, ( k„r)cos(fa) (22) 

n=l 

where Io(k n r) and I\(k n r) are modified Bessel functions, 

(23) 


where the constant B 0 (BO) is a program input. 

Figure 2 shows the geometry and axial magnetic field for 
single-period PPM focusing. The gap length between magnets 
is 2 g (TWOGCM), the inner radius of the magnet stack is a 
(ABARCM), the magnetic period is P (LMAGCM), and the 

B z (a,z) 



Ring 

magnet 



with P the magnetic period, and 



8 B 0 [fi sinh fig cos k„g + k n cosh fig sin k„g \ 

p(ji 2 + ^/o(M) 

n odd 

0, n even 


(24) 


Alternatively to entering B z (r,z) at r = a, measured on-axis 
values 5 z (0,z) over a magnetic period can be entered (as the 
input parameter BZDATA). 


B z (a,z) 



P 



^ P ► 

P ► 

Figure 2.— Geometry and axial magnetic field for single-period PPM focusing. Figure 3.— Geometry and axial magnetic field for double-period PPM focusing. 
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Figure 3 shows the geometry and axial magnetic field for 
double-period PPM focusing. The magnetic field is given by 
equations (21) and (22) with 


(l cos ^(8 B 0 ) ( fi sinh pg cos k n g + k n cosh pg sin fc*#) 


a n = 


P(ji 2 + kl)l 0 (k n a) 


n odd 


10 , 


n even 


( 25 ) 


As before, the measured on-axis values 2? z (0,z) over a 
magnetic period may be entered instead. 

Descriptions of the program input parameters are given as 
follows: 


FREQGH 

I0BMA 

JSCF 

KAV13 


Name 


Description 


ABARCM 

ACM 

ALPHL(K) 


ALPHLR(K) 


B0(K) 


BILDP(K) 

BCM 

BZDATA(I) 


DRDZ(I) 


Inner radius, a, of magnet stack from iron 
pole piece to tube axis, cm. ABARCM 
needs to be loaded only if NBZDAT = 0 
Tunnel radius, a, cm 
Power attenuation of forward wave in going 
from cavity K to cavity K + 1, dB/cavity; 
K = 1, LASTCV. To simulate a sever, a 
large value (e.g., 350 dB) is chosen for 
ALPHL in the sever cavity 
Power attenuation of backward wave in 
going from cavity K + 1 to cavity K, 
dB/cavity; K = 1, LASTCV. ALPHLR 
needs to be loaded only if KLOSS = 1 
Axial magnetic field at inner radius of 
magnet stack and midpoint of gap in K th 
magnetic section, B 0 , T; K = 1, 
LASTMG. B0(K) needs to be loaded only 
if NBZDAT = 0 

Phase shift of voltage for cavity K divided 
by 7T, (0 x L) K /r 9 rad; K = 1, LASTCV 
Initial beam radius, b 0 , cm 
Experimentally obtained value of axial 
magnetic field on the axis for NBZDAT A 
equally spaced data points. 1=1 corre- 
sponds to the beginning of a magnetic sec- 
tion, and I = NBZDATA + 1 corre- 
sponds to the end of a section. 
BZDATA(I) needs to be loaded only if 
NBZDAT *0 

Initial value of dr/dz for I th ring; 1 = 1, 
NRINGS. I = 1 corresponds to the inner- 
most ring 


KIMP 

KLMAG 

KLOSS 

KPPM 

KPRINT 

KREL 

KSMSIG 

KSPACE 


Frequency, GHz 

Beam current at tube entrance, / 0 , mA 
Number of integration steps between cal- 
culations of space-charge forces 
Number of the first cavity in the two- 
dimensional region 
KAV13 = 1: two-dimensional model 
KAV13 = LASTCV + 1: one-dimensional 
model 

KIMP = 0: Pierce impedance is entered as 
input in ZIMP 

KIMP = 1 : total impedance is entered as 
input in ZIMP 

KLMAG = 0: cavity lengths and magnetic 
sections coincide 

KLMAG = 1: cavity lengths and magnetic 
sections do not coincide 
KLOSS = 0: ALPHLR(K) is set equal to 
ALPHL(K); thus ALPHLR(K) does not 
need to be loaded. This is the usual case 
KLOSS = 1: the backward wave attenuation 
is given by ALPHLR(K) 

KPPM = 1 : single-period PPM focusing is 
to be used 

KPPM = 2: double-period PPM focusing is 
to be used 

KPRINT = 0: print input data 
KPRINT = 1 : do not print input data 
KREL = 0: use relativistic equations of 
motion 

KREL = 1 : use nonrelativistic equations of 
motion 

KSMSIG = 0: do not print small-signal 
parameters 

KSMSIG = 1: print small-signal param- 
eters 

KSPACE = 0: calculate space-charge 

forces 

KSPACE = 1 : set space-charge forces equal 
to zero 
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KSOLEN 

KSOLEN = 0: use PPM focusing 
KSOLEN = 1: use solenoid focusing; 



B z = B0, B r = 0 

NBWM 

KVEL(I) 

Number of cavity for I th printout of 



normalized axial velocities of rings; 
1=1, NVEL 

NBZDAT 

KWRITM 

KWRITM = 0: do not print normalized 
masses of rings 

KWRITM = 1: print normalized masses of 



rings whenever normalized axial velo- 
cities of rings are printed 


KWRITV 

KWRITV = 0: do not print normalized 
radial velocities of rings 
KWRITV = 1 : print normalized (to initial 

NCAVSS 


axial velocity) radial velocities of rings 
whenever normalized axial velocities of 

NDISKS 


rings are printed 

NMAX 

LASTCV 

Number of last cavity (LASTCV < 200) 


LASTMG 

Number of magnetic sections (LASTMG 



< 200). A magnetic section is one-half 
of the magnetic period for single-period 
PPM focusing and one-fourth of the 
magnetic period for double-period PPM 
focusing 

NPGRID 

LCIRCM(K) 

Length of X th cavity, L K , cm; K=l, 



LASTCV 

NPSC 

LGAPCM(K) 

Length of gap in X th cavity, 2t K , cm; 



K= 1, LASTCV 


LMAGCM(K) 

Length of X th magnetic section, cm; 



K = 1, LASTMG. LMAGCM(K) needs 

NRINGS 


to be loaded only if KLMAG = 1 

NVEL 

MSHAPE(K) 

Electric field shape parameter for X th 
cavity, m, m -1 ; K = 1, LASTCV; with 
ratio 

NXGRD1 


apedjse) = cosh( ^ } (26 ) 

E z (d,Z m center) 

NXGRD3 


where m = MSHAPE(K) and t is half- 
length in meters of X th cavity gap. (This 

NXMAG 


formulation is discussed in detail in 
reference 4.) For blunt tunnel tips, 
MSHAPE = 1 


MU(K) 

Magnetic field shape parameter for X th 
magnetic section, /x, m -1 ; K— 1, 
LASTMG; with ratio 

^(^Zgapedge) =cosh(p ^ (2?) 

B z (d,Z g a p center) 

PCDPA 


PINDBM 

TOLDV 


where jx = MU(K) and g is half-length in 
meters of X th magnetic gap. MU(K) 



needs to be loaded only if NBZDAT = 0. 
For a ratio equal to one, set MU = 1 
NBWM = 0: no backward wave 
NBWM > 0: number of last cavity con- 
sidered in calculating backward wave 
NBZDAT = 0: magnetic fields are deter- 
mined by input data BO and MU 
NBZDAT > 0: magnetic fields are deter- 
mined from experimental data. The num- 
ber of experimental data values for 
the axial magnetic field on the axis is 
NBZDAT + 1 (NBZDAT < 100) 
Number of cavity for which small-signal 
parameters are calculated 
Number of disks in beam wavelength 
(NDISKS < 24) 

Number of terms in summations stored in 
tables, excluding space-charge force 
tables and magnetic field tables 
(NMAX < 40) 

Number of grid points in radial direction 
is NPGRID + 1 . Refers only to the grid 
for calculating forces on rings due to 
electric and magnetic fields (NPGRID 
< 20 ) 

Number of grid points in radial direction 
is NPSC + 1. Refers only to grid for 
calculating forces on rings due to space- 
charge fields (NPSC < 20) 

Number of rings per disk (NRINGS < 4) 
Number of cavities for which ring axial 
velocities are printed out 
Number of grid points in axial direction 
per cavity in one-dimensional region 
(NXGRD1 < 64). NXGRD1 must be a 
multiple of four 

Number of grid points in axial direction 
per cavity in two-dimensional region 
(NXGRD3 < 64). NXGRD3 must be a 
multiple of four 

Number of grid points in axial direction 
per magnetic section (NXMAG < 64, 
normally the same as NXGRD3). 
NXMAG must be a multiple of four. 
NXMAG needs to be loaded only if 
KLMAG * 0 

Cathode flux parameter is the ratio of the 
magnetic flux linking a disk defined by 
the beam’s radius at the cathode to that 
linking a disk defined by its radius at the 
input plane of the interaction region 
Input power, P in , dBm 
Error criterion for determining whether an 
additional pass through a cavity is required. 
If 
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> TOLDV (28) 

r f, new 


TOLSC 

TOLTBL 

TWOGCM(K) 

VOB 

VJUMP(K) 

ZIMP(K) 


then an additional pass is made. Typical 
value, 0.025 

Error criterion for determining whether the 
frequency of calculating space-charge 
forces should be doubled (ref. 3). Typical 
value, 0.025 

Error criterion used in electric and mag- 
netic field tables for determining the 
number of terms to include in series 
(ref. 3). Typical value, 0.001 

Length of gap in X th magnetic section, 
2 g, cm; K = 1, LASTMG 

Beam voltage, V 0 , V 

The dc voltage jump for X th cavity, V; 
K = 1, LASTCV 

KIMP = 0: ZIMP(K) is Pierce, or beam- 
coupling, impedance in ohms for the X th 
cavity (ref. 5, eq. (10.1.19)); K = 1, 
LASTCV 

KIMP = 1: ZIMP(K) is the total inter- 
action impedance in ohms for the X th 
cavity; K = 1, LASTCV. The total inter- 
action impedance for the A^ th cavity is 
defined by 


Z K — 



(29) 


where V K is the amplitude of the voltage 
across the gap and Pf K is the average 
power flow in the forward direction 


KP Pierce impedance, fi 

ZC Total impedance, Q 

C Pierce’s gain parameter, C 

B Pierce’s fr, (u Q -v p )/v p C 

D Pierce’s loss parameter, d 

DGAIN Small-signal gain per cavity, dB/cavity 

QC Pierce’s space-charge parameter, QC/C 

A1A2 Launching loss, A x + A 2 , dB 


The following data are printed for each cavity: 


Name 


Description 


CAV 

VMAG 

ISMAG 

ISPHA 

LSGAIN 

POUT 

AVERHO 

RMSANG 


RMSVEL 


PKE 

INTRC 

PRF 

PBW 

PLC 

PBAL 


Cavity number 
Magnitude of gap voltage, V 
Magnitude of induced current normalized 
to input current 

Phase of normalized induced current 
divided by 7 r, rad 
Large-signal gain, dB 
Output power divided by I 0 V 0 
Average of the ring centroid radii divided 
by tunnel radius 

Root mean square (RMS) value of the 
angle that the velocity vectors of the rings 
make with the axis, deg 
Root mean square (RMS) value of the 
normalized radial velocities of the rings 
with respect to the initial axial velocity 
Change in beam kinetic power from initial 
beam kinetic power divided by I 0 V 0 
Power loss due to beam interception 
divided by / 0 k 0 

Power in forward wave divided by I 0 V 0 
Power in backward wave divided by IqVq 
Cumulative power loss (except for loss due 
to beam interception), divided by I 0 V 0 
Power balance equal to 


Model Output 

The program output includes three parts: the printing of 
small-signal parameters (if KSMSIG = 1), a cavity-by-cavity 
printing of selected data, and the printing of data on ring 
dynamics at selected cavities. 

The small-signal parameters (refs. 5 and 6) are as follows: 

Name Description 

U0 Initial beam velocity, « 0 , m/sec 

BEB Product of beam propagation constant and 

beam radius, (c c/u 0 )b 0 , where co is 
angular frequency and b 0 is initial beam 
radius 

BIB Product of wave propagation constant and 

beam radius, ( a)/v p )b 0 , where v p is wave 
phase velocity 


PKE - PJUMP + PRF + PBW + PLC + INTRC 

1 + (30) 

kVo 

where PJUMP is the cumulative power due 
to all the voltage jumps up to the current 
cavity 

SC Number of integration steps between cal- 

culations of space-charge forces 

At selected cavities, the following data are printed: 

Normalized axial velocities Axial velocities normalized 

with respect to the initial 
beam velocity 

Normalized radial velocities Radial velocities normalized 

with respect to the initial 
beam velocity. Printed only 
if KWRITV = 1 
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Normalized ring masses Ring masses normalized with 

respect to initial ring mass. 
Printed only if KWRITM 
= 1 

Ring radius divided by a Ring radii normalized with 

respect to tunnel radius a 


Modeling the CTS Tube 

To test the model, it was used to simulate the 
Communications Technology Satellite (CTS) traveling-wave 
tube (ref. 7), which is a 200-W, 12-GHz TWT with 58 cavities 
and 3 sections, with the output section containing 3 velocity 
tapers. The geometric dimensions and electrical parameters 
from 12.00 to 12.20 GHz in 50-MHz intervals for each of 
the 58 cavities are given in tables I to IV of reference 1. A 
sample input dataset and its resulting output are shown in the 
appendix. 

Figure 4 shows the experimental small-signal gain as a 
function of frequency at 0-dBm input power drive, 11.3-kV 
cathode voltage, and 71-mA beam current (ref. 1). Connolly 
and O’Malley compared this curve to that obtained with their 
one-dimensional model. Their results can be duplicated by 
setting KAV13 = 59, wherein this model reverts back to one- 
dimensional. The small-signal gain as a function of frequency 
both with and without a backward wave included is shown 
in figure 4. These results duplicate those of Connolly and 
O’Malley. The results without backward wave are simply those 



Figure 4.— Model and experimental results for small-signal gain as function 
of frequency for Communications Technology Satellite traveling-wave tube. 
(PINDBM = 0, YOB = 11 300, IOBMA = 71.) 


obtained after just one integration pass through the tube. To 
obtain the results with backward wave, multiple passes through 
the tube were performed until the solution converged. At the 
low end of the frequency band there was a significant 
difference between results with and without backward wave, 
and approximately ten passes were needed to obtain con- 
vergence. Through the rest of the frequency band the back- 
ward wave did not make a significant contribution, and only 
a few passes were needed. 

The one-dimensional model with backward wave simulates 
experimental small-signal gain as a function of frequency very 
well except for a ripple. It is believed that this ripple was 
primarily due to wave reflections between the output coupling 
cavity and the last sever. The model does not presently have 
the capability to model this reflection. 

To utilize the two-dimensional capability of the model, it 
is necessary to specify the magnetic field. A double-period 
PPM design is used in which there are four cavities per 
magnetic period. The inner radius of the iron magnetic pole 
pieces is equal to that of the tunnel radius and gives 
ABARCM = ACM = 0.0635. Figure 5 (from ref. 7) shows 
that the measured magnetic field at the axis was approximately 
0.15 T over the first 29 cavities and 0.22 T over the second 
29 cavities. For double-period PPM focusing, the program 
requires that the magnetic field input data be grouped in 
sections of four. Thus to model this magnetic focusing, the 
input parameter B0 (which is specified at the inner radius of 
the magnet stack) was chosen to give a maximum magnetic 
field at the axis equal to 0.1500 T over the first 28 cavities 
and equal to 0.2200 T over cavities 29 through 58. To 
determine the proper value of B0, the ratio between the 
maximum magnetic field at the magnet stack radius to that 
at the axis needs to be calculated from 


.30 



.30 1 

6 12 18 24 30 36 42 48 54 

Cavity number 


Figure 5.— Experimental axial magnetic field at axis of Communications 
Technology Satellite traveling-wave tube (from ref. 6). 
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£ a Jo (M) 

J3 0 (r = a,z = gap center) _ n=l 

B z (r = 0,z = gap center) » 

£ a » 

« = 1 

where is determined from equation (24). For this tube the 
ratio is 1.25, and thus the input is determined to be 0.1875 
T for the first 28 sections and 0.2750 T for the next 30 sections 
(BO = 28*0.1875, 30*0.2750). 

With the above magnetic field specifications, the two- 
dimensional model was run at small signal (PINDBM = 0) 
with both one ring and three rings (fig. 4). Results for both 
one and three rings were virtually identical to that for the one- 
dimensional model. This gives confidence that the model code 
is error free and also indicates that the one-dimensional model 
is adequate to determine the small-signal gain. 

Next, results at saturation at the center frequency of 12. 10 
GHz were compared. To determine the input power needed 
for saturation, runs were made scanning power input 
(PINDBM) with the one-dimensional specification. The 
resulting efficiency as a function of input power plot of figure 
6 indicates that saturation occurs for an input power of 24 dBm. 
This compares well with the experimentally obtained value 
of 23 dBm. 

At saturation (PINDBM = 24), the two-dimensional model 
with one, two, three, and four rings was compared with the 
one-dimensional model and experimental data, and the results 
are shown in table I. The results from all the models compare 
well with experiment in gain and beam efficiency. The 
experimental results indicate an interception of 6 percent, 
whereas the two-dimensional models have only 1 -percent 
interception. This is probably because thermal electrons were 
not taken into account. 

TABLE I. -GAIN, BEAM EFFICIENCY, AND 
INTERCEPTION FOR ONE- AND 
TWO-DIMENSIONAL MODELS 
AND EXPERIMENT 

[CTS TWT at saturation at 12.10 GHz.] 


Model 

Gain, 

dB 

Beam 

efficiency, 

percent 

Interception, 

percent 

One-dimensional 

29.4 

27.3 

0 

Two-dimensional 
One ring 

28.8 

23.6 

1 

Two rings 

28.8 

24.0 

1 

Three rings 

28.8 

23.8 

1 

Four rings 

29.0 

24.9 

1 

Experiment 

30.4 

24.9 

6 



0 4 8 12 16 20 24 28 32 

Input power, P in , dBm 


Figure 6.— Beam efficiency as function of input power for one-dimensional 
model with backward wave at frequency 12.10 GHz for Communications 
Technology Satellite traveling-wave tube. 

A small decrease occurs in both the gain and efficiency in 
going from the one- to two-dimensional model. This is 
primarily due to the beam interception in the two-dimensional 
model. There is virtually no difference among the results with 
one, two, and three rings. This gives added confidence that ’ 

the code is error free and indicates that the one-ring model ] 

is adequate for modeling the two-dimensional electron i 

trajectories. The accuracy is only minimally improved by \ 

increasing the number of rings. 4 

Conclusions 

The NASA large-signal, two-dimensional computer model 
for a coupled-cavity traveling- wave tube was described, and 
details of specifying the input were given. The model was 
tested by comparing its results to those of the one-dimensional 
model and experimental data for the CTS TWT. The 
agreement was excellent at both small-signal and saturation. 

It was found that the one-ring, two-dimensional model gave 
very good results, and accuracy was not improved measurably 
when more rings were used. Although the one-dimensional 
model is adequate for calculating gain and efficiency, the two- 
dimensional model computes the trajectories necessary for 
designing the magnetic focusing and calculating the beam 
entrance conditions required for multistage collector design. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, October 28, 1986 
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Appendix -Sample Dataset for CTS TWT at 12.10 GHz 


Input 

This input dataset corresponds to the two-dimensional, one- 
ring run at 12. 10 GHz of figure 4. The number of passes is 
controlled by the number of “&1NPUT &END” lines at the 


S< INPUT 

ABARCM=. 0635, 

ACM=. 0635, 

ALPHL=2*0. 1,7*. 06, 2*. 1,350. ,2*. 1, 11*. 06, 2*. 1,350. ,2*. 1, 

9*. 06, 4*. 065, 2*. 07, .075, 5*. 08, .09, . 1,2*. 11,3*. 1, 

80=28* . 1875,30*. 275, 

B1LDF‘=1 . 51 , 1.36,7*1.23, 1.36,3*1.51, 1.36, 11*1.23, 1.36,3*1.51, 
1.36,9*1.23,3*1.25,5*1.25,9*1.24, 1.36, 1.51, 

BCM=. 03175, 

DRD 2=4*0. 0, 

FREQGH=12. 10, 

I0BMA=71 .0, 

JSCF=1, 

KAV 1 3= 1 , 

KIMP-1, 

KLM A0=0 , 

KL0SS=0, 

KPPM=2, 

KPRINT=0, 

KREL=0, 

KSMSIG=0, 

KSPACE=0, 

KS0LEN=0, 

KVEL=58, 

KWR I TM= 1 , 

KWRITV=1 , 

LASTCV=58, 

LASTMG=58, 

LCIRCM=39*0. 3175, 3*0. 3048, 0. 2985, 0.2921 , 0.2794,0. 2667, 

5*0 . 2540 ,0.2413,0. 2286 , 5*0 .2159, 

LG AF‘CM=39*0 .0965,3*0.0914,0. 0899 , 0 . 0879 , 0 . 0843 , 0 . 0808 , 

5*0 . 0772 , 0 . 0739 , 0 . 070 1 , 5*0 . 0660 , 

MSHAPE=58*1 . , 


end. A negative value of ACM causes the program to 
terminate. As written, this file makes five passes through the 
tube. Changing KAV 13 to 59 results in one-dimensional 
calculations. 


MU=58*1. , 

NBWM=58, 

NBZ DAT-0, 

NCAVSS=1, 

NDISKS-24, 

NhAX=40, 

NF'GR I D=20 , 

NPSC-20, 

NRINGS=1 , 

NVEL-1, 

NXGRD1=16, 

NXGRD3=16, 

NXMAG=16, 

PCDPA-O. , 

PINDBM-O. , 

TOLD V=. 025, 

T0LSC-.O25, 

TOLTBL-1 . E— 03, 

TW0GCM=39*0. 0965, 3*0. 0914,0.0899, 0. 0879, 0. 0843, 0.0808, 

5*0 . 0772 , 0 . 0739 , 0 . 070 1 , 5*0 . 0660 , 

VOB-1 1300. , 

VJUMF-58*0. , 

Z IMP-400. ,780. , 1375. ,6*1725. , 1375. , 

780. . 2*400. , 780. , 1375. , 10*1725. , 1375. , 780. , 2*400. , 780. , 

1375. . 8*1725. , 3*1795. , 1710. , 1620. , 1530. , 1440. ,5*1350. , 

1240. . 1130.. 2*1015. , 810. , 500. , 300. &END 

?< INPUT &END 
INPUT &END 

*< INPUT &END 
INPUT 8<END 

«< INPUT &END 
INPUT ACM— 1000. 8<END 
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Output After Fifth Pass 


ORJGKU PAGE IS 

OF POOR 


CAV VMAG ISMAG TSPHA LS 

1 0.82220E+00 0.0000 -0.8377 -0 

2 0.13366E+01 0.0001 0.1510 -0 

3 0.15219E+01 0.0003 0.8426 -0 

4 0.17207Ef01 0.0007 -0.4550 -0 

5 0.19282E+01 0.0012 0.2624 -0 

6 0.19443Ef01 0.0019 0.9943 0 

7 0 . 1 7 384E* 0 1 0.0029 -0.2629 0 

8 0.17489E+01 0.0040 0.4863 0 

9 0 . 20 751E+0 1 0.0052 -0.76 38 0 

10 0.19957Ef01 0.0066 -0.0136 1 

11 0.165 1 8E+0 1 0.0081 0.7353 2 

12 0.14358E+01 0.0096 -0.5157 3 

13 0.43992E+00 0.0111 0.2332-15 

14 0.74020E+00 0.0124 0.9827 -7 

15 0. 11171 E* 01 0.0137 -0.2660 -2 

16 0.29457E+01 0.0147 0.4852 1 

17 0.44361E+01 0.0153 -0.7646 4 

18 0 . 4 9053E+ 0 1 0.0157 -0.0237 6 

2 9 0 , 4 924SE + 0 1 0.0160 0.7105 8 

20 0.61483E+01 0.0163 -0.5639 10 

21 0.79287E+01 0.0167 0.1524 11 

22 0.85851E+01 0.0174 0.8552 12 

23 0 . 8 36 68E+ 0 1 0.01 91 -0.4469 13 

24 0 . 9 2377E*01 0.0218 0.2524 14, 

25 0 . 11131 E+ 02 0.0256 0.9587 15. 

26 0.10635EF02 0.0306 -0.3295 15 

27 0 85624E+01 0.0366 0.3916 16 

28 0.72345E+01 0.0432 -0.8805 17 

29 0 . 91412E+01 0.0494 -0.1487 "2 

30 0 . 14036E+02 0.0544 0.5874 5 

31 0 . 1 0259E+02 0.0598 -0.6775 10 

32 0.22153E+02 0.0647 0.0699 14 

33 0 . 36047E*02 0.0660 0.8184 17 

34 0 . 33862E+02 0.0650 -0.4505 19 

35 0.11928E+02 0.0660 0.2654 21 

36 0.22912E+02 0.0686 0.9831 22 

37 0.48434E+02 0.0696 -0.3062 23 

38 0 . 4 9826E+ 02 0.0718 0.3832 24. 

39 0.25493E402 0.0809 -0.9302 25, 

40 0 . 30 007E+02 0.0948 -0.1995 26 

4 1 0 . 64 522E+ 02 0.1 084 0.5632 27 

42 0 , 68061 E-f 02 0.1249 -0.6787 28 

43 0.44037E+02 0.1472 0.0948 30 

44 0 . 6220 2E + 02 0.1 723 0.9014 31 

45 0.9781 3E+02 0.1954 -0.2530 32 

46 0 . 1 0238E+03 0.2146 0.6382 33 

47 0 .824 03E + 02 0.2301 -0.4240 34 

48 0.80819E+02 0.2419 0.5402 34 

49 0.96046E+02 0.2484 -0.4895 33 

50 0 . 8 97 1 3E t 02 0.2487 0.4895 31 

51 0.49928E*02 0.2427 -0.5235 30 

52 Q.15Q8SE+02 0,2319 0.4932 29 

53 0 . 58 935E + 02 0.2197 -0.4409 30 

54 0 . 7 5847E* 02 0.2074 0.6686 31 

55 0 . 67843E + 02 0.1 944 -0.2064 31 

56 0 .41894E+02 0,1788 0.9095 30 

57 0.27019E+02 0.1624 0.0196 30 

58 0 . 25125E + 02 0.14 93 -0.8767 30 

NORMALIZED AXIAL VELOCITIES 


GAIN POUT AVERHO RMSANG RMSVEL 

.0011 0.0000 0.3485 0.0063 0.0020 

.0920 0,0000 0.3345 0.0112 0.0036 

.1409 0.0000 0.3148 0.0132 0.0042 

.1109 0.0000 0.2943 0.0122 0.0039 

.0537 0,0000 0.2775 0.0085 0.0027 

.0351 0.0000 0.2690 0.0020 0.0006 

1858 0.0000 0.2714 0.0051 0.0016 

4528 0.0000 0.2841 0.0106 0.0034 

9003 0.0000 0.3031 0.0131 0.0042 

4966 0.0000 0.3239 0.0128 0.0041 

3989 0.0000 0.3418 0.0094 0.0030 

2458 0.0000 0.3523 0,0037 0.0012 

0619 0.0000 0.3530 0.0028 0.0009 

6348 0.0000 0.3439 0.0087 0.0028 

2815 0.0000 0.3269 0.0125 0.0040 

8465 0.0000 0.3062 0.0132 0.0042 

.7670 0.0000 0.2868 0.0110 0.0035 

9891 0.0000 0.2730 0.0060 0.0019 

7681 0,0000 0.2690 0.0011 0.0003 

2349 0.0000 0.2761 0.0077 0.0025 

4695 0.0000 0.2920 0.0119 0.0038 

5173 0.0000 0.3122 0.0133 0.0042 

4437 0.0000 0.3323 0.0117 0.0037 

2926 0.0000 0.3475 0.0071 0.0023 

1155 0.0000 0.3540 0.0011 0.0004 

8512 0.0000 0.3505 0.0054 0.0017 

6654 0.0001 0.3378 0.0105 0.0033 

3924 0.0001 0.3188 0.0130 0.0041 

1107 0.0000 0.2781 0.0376 0.0119 

3045 0.0000 0.2099 0.0464 0.0148 

6878 0.0000 0.1459 0.0313 0.0100 

8043 0.0000 0.1249 0.0063 0.0020 

6534 0.0001 0.1629 0.0388 0.0123 

7511 0.0001 0.2323 0.0459 0.0146 

4240 0.0002 0.2947 0.0311 0.0099 

8286 0.0002 0.3229 0.0041 0.0013 

9874 0.0003 0.3063 0.0247 0.0079 

9305 0.0004 0.2509 0.0438 0.0139 

7540 0.0005 0.1800 0.0430 0.0137 

.6363 0.0006 0.1311 0.0180 0.0057 

.6729 0.0007 0.1330 0.0205 0.0065 

.8200 0.0009 0.1847 0.0442 0.0141 

.0640 0.0013 0.2526 0.0444 0.0141 

.3741 0.0017 0.3059 0.0275 0.0088 

.6210 0.0023 0.3299 0,0077 0.0025 

.6258 0.0029 0.3241 0.0159 0.0051 

.1664 0.0033 0.2937 0.0324 0.0103 

.0884 0.0032 0.2465 0.0412 0.0131 

.3551 0.0027 0.1953 0.0381 0.0121 

.9549 0.0020 0.1563 0.0216 0.0068 

.2004 0.0013 0.1453 0.0057 0.0018 

. 3929 0.0011 0.1650 0.0275 0^0087 

.3882 0.0014 0.2053 0.0414 0.0132 

.5278 0.0018 0.2523 0.0447 0.0142 

.6560 0.0018 0.2981 0.0393 0.0125 

.8494 0.0015 0.3340 0.0274 0.0087 

.0280 0.0013 0.3548 0.0116 0.0037 

.0613 0.0013 0.3568 0.0078 0.0025 


PKE 

INTRC 

PRF 

FBW 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0001 

0.0000 

0.0000 

0.0000 

-0.0002 

0.0 0 00 

0.0000 

0.0000 

-0 . 0003 

0 .0000 

0.0000 

0.0000 

-0.0005 

0.00 00 

0.0000 

0 . 0000 

-0 .00 0 5 

0.0000 

0.0000 

0.0000 

-0 .0 0 0 5 

0.0000 

0 .0000 

0 . 0000 

-0.0004 

0.0000 

0.0000 

0.0000 

-0 .000 3 

0.0000 

0.0000 

0 .0000 

-0.0002 

0.0000 

0.0000 

0.0000 

-0.0001 

0 .0000 

0 .0000 

0 .0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 .0000 

0.0000 

0.0000 

0.0000 

-0.0001 

0.0000 

0.0000 

0 .0000 

-0.0002 

0.0000 

0.0000 

0.0000 

-0.0003 

0.0000 

0.0000 

0.0000 

-0.0004 

0.0000 

0 .0000 

0.0000 

-0.0005 

0.0000 

0.0000 

0.0000 

-0.0005 

0.0000 

0.0000 

0.0000 

-0.0005 

0.0000 

0.0000 

0.0000 

-0.0004 

0.0000 

0.0000 

0.0000 

-0.0003 

0.0000 

0.0000 

0.0000 

-0.0001 

0.0000 

0.0000 

0.0000 

-0.0001 

0.0000 

0.0000 

0 .0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0001 

0.0000 

0.0000 

0.0000 

-0.0002 

0.0000 

0.0001 

0.0000 

-0.0003 

0.0000 

0.0001 

-0 .0001 

-0.0005 

0 .0000 

0.0000 

-0.0001 

-0.0011 

0.0000 

0 .0000 

-0 .0001 

-0 .0017 

0.0000 

0.0000 

-0.0002 

-0.0020 

0.0000 

0 . 0000 

-0 .0002 

-0 .0016 

0.0000 

0 .0001 

-0 .0001 

-0 . 0010 

0.0000 

0 . 0001 

-0.0001 

-0.0006 

0.0000 

0.0002 

-0 .0002 

-0 . 0005 

0.0000 

0 .0002 

-0 .0002 

-0 .0007 

0.0000 

0 . 0003 

-0.0001 

-0.0012 

0.0000 

0.0004 

-0 .0001 

-0.0018 

0 .0000 

0 .0005 

-0.0002 

-0.0025 

0.0000 

0.0006 

-0.0002 

-0.0027 

0.0000 

0.0007 

-0.0001 

-0 . 0024 

0.0000 

0. 0009 

-0.0001 

-0.0020 

0 . 0000 

0.0013 

-0 .0002 

-0.0020 

0 .0000 

0.0017 

-0 .0003 

-0.0026 

0 . 0000 

0.0023 

-0 .0002 

-0.0035 

0.0000 

0.0029 

-0 .0001 

-0 .0041 

0.0000 

0.0033 

0 .0000 

-0 .0042 

0.0000 

0.0032 

-0.0002 

-0.0040 

0 . 0000 

0 .0027 

-0.0004 

-0 . 0035 

0.0000 

0 .0020 

-0.0006 

-0.0030 

0.0000 

0.0013 

-0 .0006 

-0.0028 

0.0000 

0 .0011 

-0.0005 

-0 . 0028 

0 .0000 

0.0014 

-0 .0003 

-0.0030 

0 .0000 

0.0018 

-0 . 0002 

-0 .0029 

0 .0000 

0.0018 

-0 .0001 

-0 . 0025 

0.0000 

0.0015 

0.0000 

-0.0022 

0 .0000 

0.0013 

0.0000 

-0 .0022 

0 .0000 

0.0013 

0.0000 


PLC PBAL SC 
0 .0000 1.000 1 
0.0000 1.000 1 
0.00001.000 1 
0.0000 1,000 1 
0.0000 1.000 1 
0 . 000 0 0 . 999 1 

0.0000 1.000 1 
0 . 0000 1.000 1 
0 .0000 1.000 1 
0 .0000 1.000 1 
0 . 0000 1 . 000 1 
0 .0000 l . 000 I 
0 .0000 1 . 000 1 
0.0000 1.000 1 
0 . 0000 1 . 000 1 
0 .0000 1.000 1 
0.0000 1.000 1 
0.0000 1.000 1 
0.0000 0.999 1 

0.0000 1.000 1 
0.0000 1.000 1 
0 .0000 1.000 1 
0.0000 1.000 1 
0.0000 1.000 1 
0.0000 1.000 1 
0.0000 1.000 1 
0.0001 1.000 1 
0.0001 1.000 1 
0 .0002 1 . 000 1 
0 .0002 0.999 1 

0.0002 0.998 1 

0.0002 0.998 1 

0.0002 0.999 1 

0 .0002 0.999 1 

0.0002 1.000 1 
0.0002 1.000 1 
0.0002 1.000 1 
0.0002 0.999 1 

0 .0003 0.999 1 

0.0003 0.998 1 

0.0003 0.998 1 

0 .0003 0.999 1 

0 .00 03 0.999 1 

0.0003 1.000 1 

0.0004 1.000 1 

0 .0004 1 . 000 1 

0 .0004 1 . 000 1 

0 .0005 0 . 999 1 

0.0006 0.999 1 

0.0806 0.998 1 

0.0007 0.998 1 

0.0007 0.999 l 
0 .0007 0.999 1 

0.0008 0.999 1 

0 .0008 1.000 1 
0 .0009 1.00 0 1 

0.0009 1.000 1 

0.0009 1.000 1 


1 .01165 

1 .01154 

1 .01053 

1.00867 

1.00639 

1.00385 

1.00111 

0 . 99819 

0.99521 

0.99240 

0.98994 

0.98797 

0.98646 

0 . 98579 

0.98621 

0 . 98772 

0.98994 

0.99295 

0.99635 

1 .00010 

1 .00369 

1.00671 

1.00916 

1 .01076 

NORMALIZED 

-0.00281 

RADIAL VELOCITIES 
-0.00301 -0,00303 

-0 . 00307 

-0 .00298 

-0.00294 

-0.00328 

-0 .00370 

-0.00390 

-0.00375 

-0.00335 

-0.00271 

-0 . 00191 

-0.00124 - 

0 .00064 

-0.00033 

-0 .00030 

-0.00038 

-0.00064 

-0.00113 

-0.00158 

-0.00189 

-0.00211 

-0.00242 

NORMALIZED 
1 . 00000 

PING MASSES 
1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 . 00000 

1 . 00 0 0 0 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1 . ooooo 

1.00000 

1 . 00000 

1.00000 

1 .00000 

1.00000 

1 . 00000 

1 .00000 

RING RADIUS DIVIDED BY 
0.33865 0.34001 

A 

0.34098 

0.34370 

0,34779 

0.34841 

0.35167 

0.35834 

0.36430 

0.36950 

0.37477 

0.37988 

0 . 38115 

0 . 38544 

0.38178 

0 . 37415 

0.36707 

0.36024 

0.35260 

0.34568 

0 . 34132 

0.33889 

0.33856 

0 .33864 
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